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Abstract 

By calculating the 0(a s ) corrections to inclusive heavy-to- light sum rules we find 
model independent upper and lower bounds on form factors for B^ttIV and B^pVu. 
We use the bounds to rule out model predictions. Some models violate the bounds only 
for certain ranges of sum rule input parameters A ~ — m\, and Ai, or for certain 
choices of model parameters, while others obey or violate the bounds irrespective of 
the inputs. We discuss the reliability and convergence of the bounds, point out their 
utility for extracting V u b, and derive from them a new form factor scaling relation. 
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1 Introduction 



A precise extraction of the Cabibo-Kobayashi-Maskawa parameters is predicated on our 
ability to control the theoretical uncertainties involved in the predictions of heavy hadron 
decays. As such, it is imperative that we perform our calculations within a systematic 
expansion. Heavy quark effective theory leads to one such approximation scheme. It leads 
to predictions for the inclusive heavy hadron decay rates and spectra p, |3| as an expansion 
in a s and Agc^/m^,, as well as predictions for heavy to heavy exclusive rates at zero recoil. 
Thus, with our present theoretical methods we are capable of extracting V c b at the level 
of a few percent. However, prospects for measuring V u b are much bleaker. The standard 
method of extraction based on the lepton endpoint spectrum will likely be plagued by large 
theoretical uncertainties while an extraction from the hadronic spectrum |J will suffer 
from the difficulty of the measurement itself. Given that exclusive modes have already 
been measured, it would be nice if we could get a handle on the exclusive form factors 
themselves. Unfortunately, our ability to calculate in QCD is limited to inclusive quantities 
due to the fact that quark-hadron duality is in itself an inclusive concept. If we attempt to 
completely eliminate the excited states from our calculation we begin to probe the dynamics 
of confinement, of which we are effectively ignorant. However, if we do not eliminate the 
excited states, but merely note that their contribution is positive definite, then we can derive 
bounds on the form factors instead of predictions. 

This idea was originally formulated by Bjorken to derive upper bounds |J on form factors 
as a function of momentum transfer q 2 , including a well-known bound on the slope of the 
Isgur-Wise function at zero recoil. In addition, Voloshin|7]] was able to find a lower bound 
at zero recoil. A few years later, Bigi et al. || showed that these bounds can be derived from 
sum rules similar to those developed for deep inelastic scattering, and are in fact the leading 
order term in a systematic expansion in a s and A^cWm^. These sum rules were then used 
to refine bounds on heavy-to-heavy matrix elements at or near zero recoil f|, ||], [l(| [ll|, |l2j, [13 



and were applied as well to heavy-to-light form factors far from zero recoil In the latter 
reference, momentum-dependent lower bounds were also developed. 

These bounds are useful for several reasons. They may be used to fine-tune and discrim- 
inate among models used extensively for CKM extractions and Monte-Carlo simulations of 
backgrounds. Given the measured q 2 dependence of a form factor, the upper and lower 
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bounds may be used to extract a value for V u b without relying on models. Since model 
predictions for B — > irlV and B — > p/z7 decay rates vary by more than a factor of four, it is 
reasonable to ask whether comparing the total rate to the integrated bounds may yield a 
more accurate extraction of V^&. We will see that the B — > 7r/z7 sum rule can place a more 
restrictive lower bound on V u b than can be obtained by the union of models. An upper 
bound may be derived as well using recent workfl5|, |16| describing the shape of form factors 
to handle the problematic zero recoil region. 



In this letter we continue the work begun in ref. [13] by including the effect of radiative 
corrections on the form factor bounds. These effects are of leading order for the lower bounds, 
and in one case dominate the upper bound. Due to the length of the calculations we only 
present partial results here and leave the more technical aspects for a subsequent, longer 
publication, where we will present the full one loop radiative corrections to the hadronic 
tensor T^ v defined by 

T^(v-q,q 2 ) = ^jd*xe- i *%B{v)\T(j»\x)r{U))\B{v)) 
= -g» v T x + v*v p T 2 - ie» va(} v a qpT 3 + q^q u T 4 

+ (q»v» + qW)T 5 , (1) 

where J is a flavor- changing current, v M is the B meson four- velocity, and g M is the dilepton 
four-momentum. These corrections are useful in their own right since they can be used to 
calculate physical rates with arbitrary cuts, via a simple numerical integration. 



Let us briefly review the derivation of the bounds. We follow the notation used in $A 
and refer the reader to this reference for details. The idea is to equate the inclusive rate, 
which is calculable within perturbative QCD using an operator product expansion, with 
the sum over exclusive states. The sum over exclusive states is then truncated, and the 
equality is changed to an inequality. More explicitly, equating partonic and hadronic sums 
over intermediate states in Eq. (H), we have 

\(H\a- J|5)| 2 + \(X\a- J\B)\* 



AM B E H e ^ H 2M B (e + E H - E 



Xj 



Here is an arbitrary four-vector chosen to pick out the form factor of interest for the 
B — > H transition (eventually we will take H to be a pion or rho), e is defined as e = 
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Mb — Eh — v ■ q, and Eh = \J M H + q 1 . The first two terms represent the local cut in 
the complex v ■ q plane corresponding to the semi-leptonic decay process. The third term 
arises from the "distant cut," corresponding to the pair production process. The sum over 
states contains the usual phase space integration / d 3 p/(2E) for each particle, while Tl x ^ H 
is shorthand for 

E' = £ (2n)W\p x + q). (3) 

To justify the use of the OPE on the right hand side of Eq. @, we perform a contour integral 
over e with a weighting function as described in [IT]. We choose the weighting function W A (e) 



to be W = 9(A — e) []. Delta has the effect of cutting off the contribution from excited states, 
which in turn improves the bound beyond what one would get trivially from insisting that 
the one exclusive mode be less than the inclusive rate. However, as discussed above it is not 
possible to choose A to be too small, since as we eliminate the contributions from excited 
states we become more sensitive to hadronization effects, and errors due to duality violation 
become important. Thus, we must choose A small enough to maximize the utility of our 
bounds, yet large enough to preserve their validity, Aqcd « A < m&, where Aqcd is 
a typical hadronic scale. Furthermore, the OPE is an expansion in Aqcd/ Eh-, so we can 
expect it to converge only when the hadronic three momentum g 3 = \q\ is sufficiently large, 
q 3 » Aqcd- 

By appropriate choice of contour we may eliminate the contribution from the unphysical 
cut. The upper bound may then be found by noticing that the excited state contribution is 
positive definite, leading to 

\(H\a- J\B)\ 2 



<Jde W A (e) a^T°J E a\ (4) 



4Mb Eh 

A lower bound can found by using the fact that 

(E 1 -E H ) ^2'\(X\a-J\B)\ 2 W A (E x -E H )< (E x - E H )\(X\a ■ J\B)\ 2 W A (E X -E H ) , 

X+H X+H 

(5) 

where E\ is the energy of the first excited state more massive than H. The contribution 
of multi-particle states with energies less than that of the first excited resonance has been 
neglected, as they are suppressed by both phase space and large- A" c power counting. We 



1 In ref. |Tl|] it was shown that for b^c the results depend rather mildly on the precise form of the weight 
function. 
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therefore have both upper and lower bounds, 

\(H\a- J\B)\ 2 



J de W A (e) a^V > ^ ^ ^ > W A (e) a»*T°f E a» 



1 - 



(6) 

AM B E H J ' - f Ex-EhY K) 

Since Eq. (|2|) forces e = E'x — Eh, we see that the last term above may be interpreted as 
the average excitation energy of the higher states contributing to the sum rule, weighted by 
W A {e) a^Ta. As emphasized in [Q, the excitation energy begins at (9(Aq CD ) and receives 
radiative corrections of order a s A. Since these two quantities are numerically comparable, 
the lower bounds are not trustworthy without the radiative corrections, which we now in- 
clude. Furthermore, given that our bounds are only valid when q 3 » Aqcd, we see that 
the lower bound will only be useful if 

(CiAqcd + C 2 a s A) q 3 < 

{Ml -M 2 H ) ~ [ ) 

Thus, we see again there is a trade off between the utility and validity of the lower bound, 

although we are helped in this case by the fact that typical hadronic masses ~ lGeV are 

numerically larger than nonperturbative condensates such as A. The range in q% for which 

the lower bounds are both useful and reliable depends on the nonperturbative corrections 



CiAqcd presented in |L4| and the radiative correction C^ci^A presented here. We will see 
that the usefulness of the lower bound depends sensitively on the numerical value of the 
heavy quark parameter A. 

The values of the parameters A and Ai have been extracted from the shape of the lepton 
end-point spectrum in inclusive semileptonic B — > Xli>i decay [[L7|, [18|] . The values found were 
A = 0.39 ± .11 GeV and Ax = —0.19 ± 0.10 GeV 2 (the uncertainty being the la statistical 
error only). The errors are strongly correlated, and therefore we will in general evaluate our 
bounds for three values of the pair (A, Ax): 

Set A = (A = 0.28 GeV, Ax = -0.09 GeV) 
Set B = (A = 0.39 GeV, Ax = -0.19 GeV) 

Set C = (A = 0.50 GeV, Ax = -0.29 GeV) (8) 

The value of A2 = 0.12 GeV 2 is determined from the B — B* hyperfme splitting. Eventually, 
A may be extracted from measurements of semileptonic — > A c decays, because all form 
factors relevant to this decay can be expressed in terms of A and one universal function[[L9|. 



and this function is furthermore determined by a single "shape parameter" [16]. Information 
on A and Ax may also be extracted from the B^Xg'y spectrum. 
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2 Heavy- to- Light Sum Rules 



Let us begin by considering the bounds on the B^ixlv form factor /+ defined by 

(7r(p') I uYb | B(p)) = (p + j/)" /+ + (P - PT /- • (9) 
The choice a = (q 3 , 0,0,v ■ q) in Eq. @ leads to the sum rule 

5 M £ ^Ud-A/M^-^-A)^^ 



+ 



Ai 
6M B 

A 2 
4M| 



Mr 

5g 3 - 5Si - 7M fl + —(6E 1 + M B ) 

Is 

Mr 

2(q 3 -E x + M B ) + —(4^ - M B ) 

Q3 



(10) 



where X and K are both functions of A, q 3 , and A. Their analytic form is rather lengthy 
and will be presented in the longer version of this article. We choose A = 1.5 GeV and 
a s = 0.32 for all plots in this section. 

Figure la shows the upper and lower bounds on the form factor / + as a function of q 2 , 
and demonstrates the importance of various perturbative and nonperturbative contributions. 
The entire kinematic range has been shown, although the bounds can only be trusted for 
q 2 ?S 18 GeV 2 . The three thin solid lines correspond to the upper bound with A and Ai 
given by data sets A, B, and C, with set A being the most constraining, and set C the least. 
It is evident the upper bound is rather insensitive to nonperturbative corrections. Neither 
is it particularly sensitive to perturbative corrections, which alter the tree level result (not 
plotted) by no more than ±15% over the entire q 2 range. The tree level lower bounds for sets 
A and B are given by the dashed lines (a) and (b), respectively. As expected, the one- loop 
contributions are comparable in size to the tree level term: they move the dashed curves (a) 
and (b) into the thick solid curves (A) and (B), respectively. The lower bound for data set 
C is not shown because it is unphysical once perturbative corrections are included (that is, 
it becomes negative). 

To assess the convergence of the perturbative and nonperturbative expansions, we look 
at the upper and lower bounds at q 2 = 17.4 GeV 2 , corresponding to qs — 1 GeV. For data 
set A with A = 1.5 GeV and a s = 0.32, the perturbative corrections in Eq. (|i"0D take the 
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Fig. la F1 9- lb 




Figure 1: Upper and lower bounds on the B^tx form factor f + as a function of momentum 
transfer q 2 , in GeV 2 . Fig. la : The three thin solid lines correspond to the upper bound for 
data sets A, B, and C as described in the text. The lower bounds for data sets A and B, 
respectively, are shown at both tree (dashed a,b) and order a s (thick solid A,B). Fig. lb : 
Thick solid lines redisplay the upper (labeled "upper") and lower (labeled "B") bounds for 
set B and the lower bound for set A (labeled 'A"). The thin lines are model predictions as 
outlined in the text. 
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values X = —0.26 and Y = —0.44. Numerically, the upper bound is 

f 2 2 

L +< !l < 0.500 - 0.055 + 0.008 - 0.021 - 0.006 

E v Mb ~ 

0(1) + 0(A) + C(Aq CD ) + 0(a s ) + 0(a s A) , (11) 

where the lower line indicates the order in the double expansion in a s and Aqcd- Both 
expansions converge quite well at this value of q$. For the lower bound, we note that 
quantities like E\ — q% — A, which represent a mismatch between perturbative and hadronic 
endpoints, are (9(Aq CD ). We evaluate such terms exactly, and for this reason do not show 
the expansion in A/q 3 , only Ai/ffe and \2/q3- The resulting lower bound is 

fill 



E n M B 



> 0.23 - 0.01 - 0.12 



The nonperturbative corrections appear to be under control. To gauge the stability of 
the lower bound with respect to perturbative corrections, it would be very useful to have 
an estimate of the two loop contribution (e.g., the a 2 /3o term). In the meantime, it is 
encouraging to note that increasing the one-loop contribution by 25% strongly modifies the 
lower bound (A) only for q 2 % 15 GeV 2 , by shifting the value of q 2 where it drops to zero 
from q 2 « 12 GeV 2 to q 2 ~ 14 GeV 2 (curve B is modified even less). 

The above analysis suggests that both the upper and lower bounds on / + are reliable 
for q 2 % 18 GeV 2 , once the nonperturbative parameters A, Ai are given. Are they useful 
as well? In Fig. 16 we again plot the two lower bounds (A) and (B) (thick solid lines), 
as well as the upper bound for the central data set B (thick solid line labeled "upper"). 
Superimposed on these bounds are four models representing the spread available in the 
literature. Three of the models include the contribution of the B* resonance by using heavy 



meson chiral perturbation theory [20|, and therefore depend on the B*-B-n coupling g and 



the B meson decay constant fs as input parameters. We adopt values used in the original 



papers. The topmost, dashed line is the prediction of Casalbuoni et al. [21 1 with g = 0.61 
and /b = -2 GeV. With this choice of input parameters, it is clearly ruled out by our upper 
bound over the entire kinematic region (and the corresponding value of V u \, eliminated). The 
thin solid line is the prediction of the Light-Front quark model by Cheung et al.f22j, with 
g = .75 and fs = -187 GeV. It violates the upper bound for q 2 > 13 GeV 2 . The dot-dashed 
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curve, corresponding to the "constrained dispersive model" of Burdman and Kamborp3 
with g = 0.5 and f B = .15 GeV, is consistent with both upper and lower bounds over 
the kinematic range where the bounds are trustworthy. Finally, the dotted curve is the 
prediction of the ISGW nonrelativistic quark model[E3. It falls below the lower bound (A) 



at q 2 = 15 GeV 2 , so it is inconsistent with the values A = .28 GeV, Ai = -0.09 GeV 2 . It 



has been pointed out[^5| that the ISGW model fails to include the contribution from the 
B* pole at large q 2 ; if modified to include this, it may well rise above our lower bounds. 
The models plotted in Fig. lb are merely a representative range - there are many other 
models and calculations. For example, the prediction of Bagan et aZ.[2?J is similar to that of 
Burdman and Kambor for q 2 ^5 18 GeV 2 , so it is also consistent with the bounds. Lattice 
simulations^] are consistent as well. 

While the upper bounds are fairly insensitive to A, the lower bounds can be somewhat 
improved by smearing over a smaller range of excitation energies. Taking A — > 1 GeV lowers 
the value of q 2 where the lower bound vanishes by 1 to 2 GeV 2 . Conversely, increasing A— > 
2 GeV weakens the lower bound, but only about half as severely as when A is decreased. We 
retain A = 1.5 GeV as a reasonable compromise between maximizing utility and reliability. 

We turn now to the B — > plv form factors a± defined by 



(pip') I uYlsb | B(p)) = fe*v + (p + pT a + + (p- p'f a. 



p ■ e 



(13) 



Using a = (E p , 0, 0, — ^3) in Eq. 



^(E p -q 3 f 



+ Ai 

+ A 2 
> 



gives the upper bound 
E 



Ai 


<?3 


13M B 


A 2 


' Q3 






M B q 2 M 2 



p_ 

3<?3 
2<?3 



El 



3qsM J 
M B q 3 



B 



a + + ~j} a+ + a 



+ 0{a s ql) 
M B E P 



M 2 



- 1 



(14) 



which is interesting because the leading order term vanishes as ^3 becomes large. As pointed 
out in [[14]], this implies the a s contribution dominates at small q 2 . 

It is worth noting that Eq. (0) implies nontrivial scaling relations for the form factors 
a±. In the limit that Aqcd << Q3, Mb, they must fall off at least as fast as 

1 , Ma) <? 3 ] 1/2 



/2 



(a+ + a_ 



(M B q 3 y 
A 2 

ii QCD 

(M B q 3 f/ 2 



M b /2 A QC b 



+ 



a 



,(g 3 ) 1/2 A 



QCD 



(15) 



The scaling with Mb follows from considering the heavy 6-mass limit, which leads toQ a_ = 
— CL|- [1 + 0(l/m.f,)]. That a + , modulo perturbative contributions, falls off with q 3 is no 
surprise, since even a simple pole goes like 1/53. However, the q^ behavior of (a + + a_) is, to 
our knowledge, a novel result. At q 2 = 0, ^3 ~ Mb and the 6-mass scaling relation need not 
hold. Instead, the full scaling relations take over. Practically speaking, the large factor 
(M B E™ ax — M 2 ) / (2M 2 ) m 12 in Eq. (|4|) leads to very tight constraints on the combination 
a + + a_, independent of the validity of the 6-mass scaling relations. 

The upper bounds on the magnitude of the dimensionless form factor combination 



A + = {M B + M p 



1 . / MbE p 

2 (fl++a - ) h^- 1 



(16) 



are plotted in thick solid lines in Fig. 2. The curves (a), (b), and (c) correspond to data sets 
A, B, and C, respectively. The parameter Ai plays an important role in this case because the 
leading order (E p — q^) 2 term suffers from a severe cancellation. In fact, for g 3 >> M p , this 
term is formally 0(Mp/q 2 ), which is subleading to the contributions from Ai and A2 (note 
that A does not enter the upper bound at tree level). The thin solid line is the tree level 
upper bound for Ai = —.19 GeV 2 (data set B). As expected, the perturbative corrections 
dominate at small q 2 . Even so, changing the one-loop contribution by 25% alters the upper 
bounds by no more than « 10% (the square root helps), so they are reasonably stable against 
perturbative corrections. Again, a two-loop calculation or estimate would be quite welcome 
here. 

Also plotted in Fig. 2 are three predictions for A + taken from the literature. The dashed 
line is the relativistic quark model of WSBp8|, which exceeds the upper bound by roughly 
a factor of five. The reason it violates the bound so spectacularly is that it does not obey 
the scaling law a_ = — a + . Indeed, in order to comply with the bounds, the scaling violation 
at q 2 = would have to be reduced by a factor of eight, (a+ + aJ)/a + < ~. The dotted 



curve is the prediction of another relativistic quark model due to Melikhov and Nikitin[^£|. 
using their "Set 1" parameters. It is consistent with the upper bounds for A x > —0.19 GeV 2 
(their other three sets are consistent for Ai > —0.09 GeV 2 ). The nonrelativistic ISGW quark 
model j24|, |30| obeys the bounds regardless of the data set. We note that at q 2 = 0, the model 
of Cheng et al. ||31|| gives A + = 0.12, which is consistent with our bounds. 

2 This relation was used in |l4| to simplify Eq. (|l4|). This is only valid for b—>c decays, since otherwise 
the smailness of (a+ + a_) is compensated by Mb/M p . 
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Figure 2: Upper bounds on the B — > p + l~V form factor A + for data sets A, B, C are given by 
the thick solid curves (a), (b), (c), respectively. The tree-level upper bound for set B is given 
by the thin solid curve. The dashed, dotted, and dot-dashed curves are the model predictions 
of WSB, MN, and ISGW, respectively, as described in the text. 
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Figure 2 illustrates the utility of the inclusive sum rules for constraining, or even elimi- 
nating, models. The sum rules can also yield bounds on the nonperturbative parameters A 
and Ai|| H [□], [32|| . For example, the condition \A + \ 2 > in Eq. ([LI]) gives a lower bound 
on Ai comparable to f52]| , if we use A = 1.5 GeV. The lower bound on A + also yields a 
restrictive sum rule, but higher order nonperturbative terms become important (this is why 
we have not plotted the lower bound) . The role of heavy-to-light sum rules in constraining 
nonperturbative parameters will be explored in greater detail in the long version of this 
article. 

3 Conclusions 

We have presented one-loop improved upper and lower bounds on the B — > irlu form factor 
/+, as well as an upper bound on the combination of B — ► plV form factors A + defined in 
Eq. (|T6|). We used the upper bounds on /+, which are exceptionally stable with respect to 
perturbative and nonperturbative corrections, to rule out certain input parameters in two 
models. A third model may be consistent with the lower bounds on / + , depending on the 
values of the measurable quantities A and Ai. The lower bounds change little if the one-loop 
corrections are increased by 25%, but an estimation of the two-loop correction is needed to 
be certain they are reliable. 

The upper bounds on A + are moderately stable against further perturbative corrections, 
but an estimate of two-loop effects would again be most welcome. The bounds are partic- 
ularly interesting because they strongly constrain the allowed size of the scale violating (in 
rrib and q^) quantity a + + a_. Indeed, the model of Wirbel, Stech, and Bauer exceeds the 
upper bound by a factor of five. A model by Melikhov and Nikitin obeys the bounds, but 
depending on the physical value of Ai, their input parameters may be constrained. The 
ISGW model, by contrast, obeys the bounds rather handily. 

The bounds given here demonstrate how inclusive heavy-to-light sum rules can constrain 
model parameters, or in some cases even invalidate their predictions. They should be useful 
not only in discriminating between models, but in constructing, constraining, and fine-tuning 
them. A variety of bounds on various combinations of form factors can be made. The in- 
clusion of a s perturbative corrections now allows us to consider, with some reliability, lower 
bounds on form factors, upper bounds on heavy quark violating form factor combinations, 
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and upper bounds on nonperturbative condensates. The results presented here will be ex- 
tended in a forthcoming publication including analytic expressions for the perturbative 



corrections. Improvements to existing sum rules may be possible by either further restricting 
the spin-parity of contributing intermediate states, or by using data to include the contri- 
butions of higher states to the hadronic sum. They should be of use as well in extracting 
model-independent limits on the CKM element V u b. 
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